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The influence of fixed temperature and fixed heat flux thermal boundary conditions on 
rapidly rotating convection in the plane layer geometry is investigated for the case of 
stress-free mechanical boundary conditions. It is shown that whereas the leading order 
system satisfies fixed temperature boundary conditions implicitly, a double boundary 
layer structure is necessary to satisfy the fixed heat flux thermal boundary conditions. 
The boundary layers consist of a classical Ekman layer adjacent to the solid boundaries 
that adjust viscous stresses to zero, and a layer in thermal wind balance just outside the 
Ekman layers adjusts the temperature such that the fixed heat flux thermal boundary 
conditions are satisfied. The influence of these boundary layers on the interior geostroph- 
ically balanced convection is shown to be asymptotically weak, however. Upon defining 
a simple rescaling of the thermal variables, the leading order reduced system of govern¬ 
ing equations are therefore equivalent for both boundary conditions. These results imply 
that any horizontal thermal variation along the boundaries that varies on the scale of 
the convection has no leading order influence on the interior convection. 
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1. Introduction 

One of the simplest and most commonly studied systems for investigating convection 
dynamics is the so-called Rayleigh-Benard configuration, consisting of a Boussinesq fluid 
layer of depth H confined between plane-parallel boundaries, and heated from below. 
The constant gravity vector g = —gz points vertically downwards. Two limiting cases 
for thermal boundary conditions are often considered when posing the problem math¬ 
ematically: (1) ‘perfectly conducting’, or fixed temperature (FT), boundary conditions 
in which the temperature is held fixed along the bounding surfaces; and (2) ‘perfectly 
insulating’, or fixed flux (FF), boundary co nditions in whic h the norma l derivative of the 
temperature is fixed at the boundaries (e.g. Chapman fc Proctoilll98C)h . Thermal bound¬ 
ary conditions of geophysical and astrophysical relevance are often considered to reside 
somewhere between these fixed flux and fixed temperature limits. 

For a Newtonian fluid of constant thermal expansivity a, kinematic viscosity z/, ther¬ 
mal diffusivity n. the non-dimensional Rayleigh number quantifies the strength of the 
buoyancy force. For the FT and the FF cases we have, respectively, 


Rclft 


■g \ I R 


VK 


Rclff = 


agPH 4 


( 1 . 1 ) 
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where AT is the fixed temperature difference between the top and bottom boundaries and 
/3 is the fixed temperature gradient maintained at the boundaries. The Prandtl number 
quantifies the relative importance of viscous and thermal diffusion as Pr = v/ k. Upon 
defining the non-dimensional measure of heat transfer via the Nusselt number, 


Nu = 


total heat transfer /3H 
conductive heat transfer AT ’ 


( 1 . 2 ) 


it is straightforward to show that the two Rayleigh numbers defined above are related 
simply by Rapp = NuRapp- We thus see that for linear convection in which Nu = 1 
the two Rayleigh numbers are equivalent. For nonlinear convection in which the critical 
Rayleigh number has been surpassed, Nu > 1 is achieved by adjustment of the temper¬ 
ature gradient ft at fixed AT for FT boundaries, and vice versa for FF boundaries. 

Linear stability shows that for the case of non-rotating co nvection the most un - 
stable wavenumber is finite for FT boun dary condit i ons ( e.g. IChandrasekhar 1961 1. 
but is zero for FF boundary conditions ( Hurlc et al. 119671 1. Although previous work 


suggests that these d i fferen ces for linear convection also hold for nonlinear convection 
( Chapman fc Proctodll980l) . numerical simulations of convection show that the statistics 
for the tw o cases converge as the R ayleigh number is increased and the flow becomes 
turbulent ( Johnston fc Docrind 2009i l. 

When the system is rotating with rotation vector Cl = flz, the Ekman number, 
Eh = v/2iIH 2 , is an additional non-dimensional number required to specify the strength 
of viscous forces relative to the Coriolis force. The rapidly rotating, quasi-geostrophic 
convection limit is characterized by Ep — > 0. As of this writing, only two investigations 
of FF boundary cond itions for the ro tatin g plane layer g e ometr y have been published 
in the literature, with iDowlind ( 1988h and Takehiro et all ( 2002 ) exam i ning the weakly 


rotating and rapidly rotating linear cases, respectively. Takehiro et al. ( 20021 ) utilized a 


modal truncation approach to show that the critical parameters for the two cases should 
converge as Ep —> 0; the present work confirms this suggestion. 

In the present work we distinguish between ‘interior’ and ’boundary layer’ dynamics, 
and show that the interior governing equations are identical for the two different thermal 
boundary conditions upon a simple rescaling of the Rayleigh number and temperature. 
Because the Ep — > 0 limit is a singular perturbation of the governing equations, the 
interior equations cannot satisfy the FF boundary conditions at leading order; a dou¬ 
ble boundary layer structure is necessary to adjust both the horizont al viscous stresses 
and t he normal derivative of the temperature fluctuation to zero (c.f. iHeard fc Vcronis 
19711) . It is shown that the boundary layer corrections are asymptotically weak, how¬ 
ever, showing that to leading order the interior quasi-geostrophic convection dynamics 
are equivalent for both thermal boundary conditions. 

In section [2] we present the linear stability of the full Boussinesq Navier-Stokes equa¬ 
tions. In section[3]we present the asymptotic reduction of the Navier-Stokes equations in 
the rapidly rotating limit and concluding remarks are given in section |4] 


2. Linear stability of the Navier-Stokes equations 

In the present section we briefly present the linear stability of the Boussinesq Navier- 
Stokes equations for both FT and FF thermal boundary conditions. Upon scaling lengths 
with the depth of the fluid layer H and time with the viscous diffusion time H 2 /p, the 
linear system becomes 
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d t d' -w= VV, (2.2) 

Pr 

V-u = 0, (2.3) 

where the velocity vector is denoted by u = (u, v, w), and the temperature is decomposed 
into mean and fluctuating variables according to d = d + d'. For both sets of thermal 
boundary conditions d = 1 — z and the fluctuating thermal boundary conditions therefore 
become 

d' = 0, at z = 0,1, (FT) (2.4) 

d z d'= 0, at z = 0,1. (FF) (2.5) 

Stress-free, impenetrable mechanical boundary conditions on the top and bottom bound¬ 
aries are assumed throughout and given by 


w = d z u = d z v = 0, at z = 0,1. (2.6) 

The system flUD-flgp is discretized in the vertical and horizontal dimensions with 
Chebyshev polynomials and Fourier modes respectively, and formulated as a general¬ 
ized eigenvalue problem. We solve the system in primitive variable form and enforce 
boundary conditions via the tau method. The eigenvalue problem is solved with Mat- 
la b’s ‘sptarn’ fun ction. For further details of the numerical methods the reader is referred 


to lCalkins et al\ ( 20131 . 


Figure Q] shows results from the linear stability calculations. Results are given for 
both steady (Pr = 1) and oscillatory (Pr = 0.1) convection; we note oscillatory con- 
vection does not exis t for Pr ^ 1 and becomes the primary instability for Pr < 0.68 
( ChandrasekhaiHl961 ). For Eh < 10 -5 , both the asymptotically scaled critical Rayleigh 
number Fa c F 4 / 3 (Figure[TJi) and wavenumber fc c F 4 / 3 (FigureH}]) obtained from FT and 
FF boundary conditions are observed to converge to nearly equivalent values. The open 
circle shows the Ekman number Eh = 0.0745 calculated bv lDowlind (1988) at which the 
instability becomes characterized by k c ^ 0 with Ra c = 341.05. 

Figure [I]/ shows the horizontal ( x ) velocity eigenfunction for Eh = 10 -4 where Ekman 
layers can be seen at the top and bottom boundaries for the FF (dashed curve) case; a 
magnified view of the bottom Ekman layer is shown in the inset figure. The temperature 
perturbation eigenfunctions shown in FigureQJl show that both the FT and FF cases have 
identical structure in the fluid interior, whereas a thermal boundary layer is observed in 
the FF case. In the following section we present the asymptotic reduction of the Navier- 
Stokes equations to better understand and quantify this behavior. 


3. Asymptotics 


To proceed with the asymptotic development, we follow the work of ISprague et al 


( 20061) and write the governing equations using a generic non-dimensionalization such 
that 


D t u + -|-z x u = —EuS/p + T#z + — V 2 u, 

Ro Re 

(3.1) 

V • u = 0, 

(3.2) 

Dtd = P L^ 

(3.3) 


where D t (-) = dt(-) + u ■ V(-) and the velocity, pressure and temperature are denoted 
by u, p , and d , respectively. The above system has been non-dimensionalized utilizing 
the velocity scale [/, length F, time L/U 1 pressure F, and temperature T. For the FT 
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Figure 1: Linear stability of the Navier-Stokes equations for fixed temperature (FT) and 
fixed flux (FF) thermal boundary conditions, (a) Asymptotically scaled critical Rayleigh 
number and (b) critical wavenumber as a function of the inverse Ekman number for both 
steady (Pr = 1) and oscillatory (Pr = 0.1) convection, (c) Horizontal velocity and (d) 
temperature eigenfunctions for Pr = 1 and Eh = 10 -4 ; the inset figure in (c) shows a 
magnified view of the Ekman layer along the bottom boundary. In (a) and (b) the Ekman 
number of Eh = 0.0745 at which the critical wavenumber for FF b oundary conditi ons 
becomes non-zero is shown by the open circle, as first determined bv [Dowling (1988). 


and FF cases the temperature scale becomes AT and NuAT, respectively. The Rossby, 
Euler, buoyancy and Reynolds numbers are defined by 


Ro = 


U 

2QL 1 


Eu = 


P 

~PoU*' 


r = 


gaT 
U 2 ’ 


Re = 


UL 


(3.4) 


In the present work we are interested in the e = Ro —> 0 limit. In the fluid interior we 
employ multiple scales in the axial space direction and time such that 


d z ->• d z + ed z , i d t ^d t + e 2 d T , 


(3.5) 


where Z = ez is the large-scale vertical coordinate and r = e 2 t is the ‘slow’ timescale. 
It has been shown that the following distinguished limits can be taken 
governing equations to accurately model quasi-geostrophic convection (e.g. 


to reduce the 
Sp rag ue et aL 
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200 6) 


Eu=\, T = —, Re = 0(1), Pr = 0(1), 
e A e 

where T = 0(1). Scaling the velocity viscously such that U = v/L we have 

E^Ra 


= s 1 / 3 , r = 


Pr 


Re = 1, 


(3.6) 


(3.7) 


where the L-scale Ekman number is related to the H -scale Ekman number via E = 
Ejje~ 2 , i.e. L = E X ! 3 H. We keep the notation for the Rayleigh number generic in the 
sense that Ra denotes either Rapr or Roff depending upon the particular boundary 
conditions employed. Hereafter, we define the asymptotically reduced Rayleigh number 
as Ra = E 4 / 3 Ra. 

We ut il ize a composite as ymptotic expansion approach (e.g. Navfchl [20081) and, follow¬ 
ing Heard Sz Veronia (1971), decompose each variable into interior (i), middle (m), and 
Ekman layer (e) components. For instance, the dependent variable / can be written as 


/ = / W (x, Z, t, t ) + / (m) (x, y, t) + / (e) (x, y, ry, t), 


(3.8) 


where £ = z and ry = e -1 / 2 z are boundary layer variables. The above representation 
ensures that each dependent variable is uniformly valid throughout the domain. The 
boundary layer variables consist of a sum of contributions from the top and bottom 
boundary layers; for brevity, we focus on the bottom boundary layers. In the present 
work we make use of the following limits and notation 

lim (/<"*>) = (/("*)) W = Inn (>)) = ( ° = 0, (3.9) 


lim (f 1 ' 
V 


= lim (/«) = / (i) (Z = 0) = / W (0). 


(3.10) 


We then expand each variable in a power series according to 

/ W (x, Z,t) = /d ?) (x, Z,t,r) + e 1/2 f^ 2 (x, Z,t,r) + ef[ l) (x, Z,t,r) + 0(e 3/2 ). (3.11) 

Each dependent variable is further decomposed into mean and fluctuating components 
such that 

/«(x, Z, t, t) = 7 (i) (Z, r) + f'W (x, Z, t, r), (3.12) 

where the averaging operator is defined by 

1 / 


f(Z, t ) = lim 

M,[A]-)-oo [r][H] J[ t ]\a] 


f dxdy, and f = 0. 


(3.13) 


3.1. The interior equations 

By substituting decompositions for each variable of the form (13.81) into the governing 
equations and utilizing the limits (I3.9H - (13.101) . equations for each region can be derived; 
expansions of the form (13.111) are then utilized to determine the asymptotic behavior of 
each fluid region. Because the derivation of the interior equations has been given many 
times previously, we present only the salient features and direct the reader to previous 
work (e.g. Sprague et all 12006 ) for details on their derivation. The main point is that the 
interior convection is geostrophically balanced and horizontally divergence-free to leading 
order 

Z X u^ } = -V_LpW, Vx-Uq*j_ = 0 , 


( 3 . 14 ) 
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where V_l = (d x ,d y , 0). The above relations allow us to represent the geostrophic velocity 
via the geostrophic streamfunction ip^ = such that = —V x ip^z. The ver- 

tical vorticity is then Q ; ■ The interior vertical vorticity, vertical momentum, 

fluctuating heat, and mean heat equations then become 

Dt - dzw' 0 (i) = Vi Co 

yw+vi^o 




dztp® = |V (i) ■ ^ 2 

Pr 


dr 4 


,,'W 


5 (i) 


1 


= p- r ^in 


/« 




(3.15) 

(3.16) 

(3.17) 

(3.18) 


where Dp (•) = <9t(-) + u- Vj_(-) . The mean interior velocity field is zero and the mean 

momentum equation reduces to hydrostatic balance in the vertical, dzp\P = (Ra/Pr)d q ^. 

The interior system is fourth order with respect to the large-scale vertical coordinate Z. 
Two boundary conditions are supplied by impenetrability such that w'q 1 ' 1 (0) = Wq^ (1) = 
0. Although no ^-derivatives with respect to t>'P > are present in equation (13.171) . eval¬ 
uating this equation at the boundaries shows the FT conditions i?^(0) = i//*'(l) = 0 
are satisfied implicitly for the fluctuating temperature. Evaluating equation (13.161) at 
either the top or bottom boundary with the use of impenetrability shows that stress-free 
boundary conditions are implicitly satisfied as well since dzip q (0) = dztpo\ 1) = 0. 


For the case of FT thermal boundary conditions, we have 

y } (0) = 1, and y } (l) = 0. (FT) 


(3.19) 


Thus, for the FT case the boundary layer corrections are identically zero and the above 
system is complete. Numerous investigations have used the above system of equations to 
investigate rapidly rotating convection in the presence of stress-free mechanical bound¬ 
ary conditions and have shown excellent agreement with direct numerical s imul ations 
(DNS) of the Navie r-Stokes equations and laboratory experiments (|Stellmach et nl.ll2014 : 
Aurnou et al. 120151) . 

For the FF case the mean temperature boundary conditions become 


5 (0 


■(*)/ 


dz#o (0) = —1, and dzv 0 (1) = —1. (FF) 


(3.20) 


We further require dz’d'^ (0) = dzd ^(l) = 0; boundary layer corrections are therefore 
required since these conditions are not satisfied by (|3.17l) . In the following two subsections 
we determine the magnitude of these boundary layer corrections. 

3.2. The middle layer equations 

The first non-trivial fluctuating middle layer momentum equation occurs at 0(e) to yield 
the thermal wind balance 


q'OA-n — 


Z X u. 


l(m] 


= ^P't ] + 


(3.21) 


such that V_l-U 2 = 0 and w'p'"' 1 = 0. The mean velocity field u <m ) = 0. 
The leading order temperature equation for the middle layer is 




to 2 


A m ) , , '(*) 


U n 


'(0)-V_Li?2 (m) = y^-V 2 i?2 (m) , 

Pr 


(3.22) 
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with corresponding boundary conditions 

dz^ (0) + d^ m) (0) = 0, $2 (m) (£ -t oo) 0. (3.23) 

We find the first non-trivial mean temperature to be of magnitude 0(e 5 ) and therefore 
omit any further consideration of this correction. 

The first three orders of the stress-free mechanical boundary conditions along the 
bottom boundary become 

(0) = 0, dzu^l _l( 0) = 0, ^zu'/^O) + + ^Ug (e) (0) = 0. (3.24) 

Thus, the first two orders of the interior velocity satisfy stress free conditions on their 
own and therefore need no boundary layer correction. Here we have rescaled the Ekrnan 
layer velocity according to = e 5//2 Ug C - ) ; this rescaling is simply highlighting the fact 
that the Ekrnan layer velocities are significantly weaker than those in the interior. 


d z™ o,l 


3.3. The Ekrnan layer equations 


The Ekrnan lay er equations have been studied in great detail in previous work (e.g. 
Greenspan 1968ll . so we simply state the leading order continuity and momentum equa¬ 


tions as 


Vj_-Uo (e) + d v w'i e ^ =0, z x Ug e) = <9 2 Ug e) (3.25) 

2 

where w'^ = e 3 w'i e ^. All of the mean Ekrnan layer variables can be shown to be zero. 
2 

A key component in the present analysis that differs from previous work is the middle, 
thermal wind layer that enters the Ekrnan layer solution via the stress-free boundary 
conditions (13.241) . Utilizing the thermal wind relations for the middle layer that follow 
from equation (13.211) . 


= ^ r dA m \ d^ m) = §^d x 4 m) , (3.26) 

the stress-free boundary conditions along the bottom boundary can be written as 

dzu^liO) + |^V x 4 m) (0) + ^ (e) (0) = 0, (3.27) 

where = (—<9 y , d Xl 0). Solving the Ekrnan layer momentum equations for the horizon¬ 
tal components of the velocity field with the additional requirement that ) —> 0 

as £ —> oo, the continuity equation is then used to find the Ekrnan pumping velocity 


~/(e) 

Wi = 


0zd°( 0) + ^Vi4 m) (0) 


-V ( 77 

^ cos —= 

Vv^ 


(3.28) 


Thus, vertical velocities of magnitude 0(e 3 ) are induced by FF thermal boundary con¬ 
ditions and result from both finite viscous stresses within the fluid interior and hor¬ 
izontal variations of the temperature within the middle layer . This finding is closely 
analogous to the Ekrnan pumping effect first reported bv iHidd (19641) for shallow layer 
quasi-geostrophic flow in the presence of lateral temperature variations along a free sur¬ 
face. Evaluating equation (13.281) at 77 = 0 provides a parameterized boundary condition 
for the effects of Ekrnan pumping. 

The small magnitude of the Ekrnan pumping velocity (13.281) results in very weak 0(e 5 ) 
temperature fluctuations within the Ekrnan layer. Because of this, the dominant correc¬ 
tion of the FF thermal boundary conditions occurs within the middle layer and we do 
not consider the Ekrnan layer temperature any further. 
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3.4. Synthesis 

The thermal boundary layer correction given by equation (13.231) is passive in the sense 
that $2 can be calculated a posteriori with knowledge of . Thus, the leading order 
interior dynamics are insensitive to the thermal boundary conditions. The Ekman layer 
analysis shows that the first six orders of the interior vertical velocity satisfy the impene¬ 
trable mechanical boundary conditions wf l \ 0) = 0, for i = 0,..., 5/2. At 0(e 3 ) we have 
the Ekman pumping boundary conditions 

«4 W (0) = -«4 (m) (0) - dzCi\0) - ^Vi4 m) (0), (3.29) 

where we have used the Ekman pumping relation (13.281) evaluated at 77 = 0. From the 
standpoint of linear theory, the first correction to the critical Rayleigh number will there¬ 
fore occur at 0(e 3 ); this explains the linear behavior previously discussed in section [2] 


4. Discussion 

In light of the boundary layer analysis, we conclude that the leading order quasi- 
geostrophic dynamics are described by equations (I3.15I) - (I3.18I) for both FT and FF ther¬ 
mal boundary conditions. Indeed, inspection of the system shows that it is invariant 
under the following rescaling of the Rayleigh numbers and temperature variables, 


Roft = 


Roff 

Nu 


^1 .FT = Nu $ 1 ,FF 1 $0,FT = Nu $ 0,F.F • 
Integrating the time-averaged mean heat equation with respect to Z yields 

Pr (w' ( V/Jy) = dzd^pT + Nu, (FT) 


(4.1) 


(4.2) 


Pr (w' 0 (i) = d z d^ FF + 1, (FF) (4.3) 

for the FT and FF cases, respectively. The appropriate thermal boundary conditions 
have been applied at Z = 0 in the above relations. Taking either equation (14.2[> or (14.31) 
and utilizing m shows that the mean interior temperature gradient is described by 
identical equations for the two cases. This leading order correspondence is the result of 
the anisotropic spatial structure of rapidly rotating convection. 

The above results indicate that the findings of previous work on low Rossby number 
convection employing FT thermal boundary conditions can be accurately applied to the 
case of FF th ermal b oundary conditions by use of the rescalings given by equations dUD. 
Julien et ail ( 2012 fob identified four flow regimes that occur in rapidly rotating convec¬ 
tion as a function of the Prandtl and (FT) Rayleigh numbers. The so-called ‘convective 
Taylor column’ (CTC) regime is distinguished by coherent, vertically aligned convective 
structures that span the depth of the fluid. Figure [2ja) shows a volumetric rendering 
of the temperature perturbation for Pr = 7 and Roft = 46.74, or Roff = 1000 
and Nu = 21.39. The CTC regime occurs over the FT Rayleigh number range of 
20 < Rq ft 5. 55, correspo nding to a FF Rayleigh number range of 82 < RapF 1656 
(e.g. see lNieves et a/.ll2014l) . Figure Ob) shows mean temperature profiles obtained uti¬ 
lizing the FT and FF thermal boundary conditions, along with the remapped FF mean 
temperature profile. The Nusselt number Nu = 21.39 corresponds to a mean temperature 
difference of 0.0468 between the top and bottom boundaries for the FF case. 

Of particular interest in convection studies is the dependence of the heat transfer scal¬ 
ing with the strength of the thermal forcing input via Nusselt-Rayleigh number scalings 
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(a) 



(b) 


Figure 2: (a) An example volumetric rendering of the temperature perturbation from a 
simulation of the quasi-geostrophic convection equations showing the ‘convective Taylor 
column’ (CTC) regime, (b) Mean temperature profiles obtained with both FT (solid 
blue) and FF (dashed black) boundary conditions, and the rescaled FF temperature 
profile (red open circles). The parameters are Pr = 7, Roft = 46.74, Roff = 1000, and 
Nu = 21.39. 


of the form Nu ~ Ra FT . With the rescaling given in m the FF equivalent of this 

relation becomes Nu ~ Ra e p wher e (3 = a/[a + 1). For the CTC regime the expo¬ 
nent is a ~ 2.1 (jjulicn e t a/.ll2012M) . yielding /3 ss 0.68. Additionally, the final regime 
of geostr ophic turbulence ac hieves a dissipation-free scaling law with a = 3/2 such that 
/3 = 3/5 ( Julien et al. 2012 a !). Similarly, the dependence of all other variables of interest 
on the Rayleigh number (e.g. mean temperature gradient, vorticity, etc.) can also be 
remapped to the case of FF thermal boundary conditions. 


5. Conclusion 

In the present work we have shown that the leading order dynamics of rapidly rotating 
convection in a plane layer geometry are equivalent for both FT and FF thermal bound¬ 
ary conditions. FF thermal boundary conditions give rise to a double boundary layer 
structure in the limit of rapid rotation that is asymptotically weak. Our findings sug¬ 
gest that all previous work employing FT thermal boundary conditions also accurately 
describes FF thermal boundary conditi ons with i n the reg ime of asymptotic validity, i.e. 
Ro ~ A 1 / 3 <C 1 and Ra < 0{E~ 1 ^) (j.Tulien et ai.ll2012flh . This result adds to the ro¬ 
bustness of the reduced, quasi-geostrophic model in geophysical and astrophysical appli¬ 
cations where stress-free boundary conditions are typically assu med; recent wo r k ha s also 
extended the model to the case of no-slip boundary conditions ( Julien et al. 201 51) that 


are of relevance for lab oratory experiments and planetary interiors (e.g. lAurnou et al 
201 fit ICheng et flJ.11201 5l) . 


An interesting consequence of these findings is that any horizontal thermal variation 
along the boundaries that varies on the scale of the convection has no leading order 
influence on the interior convection. This finding helps to explain the results of previous 
spherical convection studies investigati ng the role of temperature variations along the 
outer boundary (e.g. Davies et al. 20091 ). 






































10 


M. A. Calkins, K. Hale, K. Julien, D. Nieves, D. Driggs and P. Marti 

REFERENCES 

Aurnou, J. M., Calkins, M. A., Cheng, J. S., Julien, K., King, E. M., Nieves, D., 
Soderlund, K. M. & Stellmach, S. 2015 Rotating convective turbulence in Earth and 
planetary cores. Phys. Earth Planet. Int. . 

Calkins, M. A., Julien, K. & Marti, P. 2013 Three-dimensional quasi-geostrophic convection 
in the rotating cylindrical annulus with steeply sloping endwalls. J. Fluid Mech. 732 , 214- 
244. 

Chandrasekhar, S. 1961 Hydrodynamic and Hydromagnetic Stability. U.K.: Oxford University 
Press. 

Chapman, C. J. & Proctor, M. R. E. 1980 Nonlinear Rayleigh-Benard convection between 
poorly conducting boundaries. J. Fluid Mech. 101 , 759-782. 

Cheng, J. S., Stellmach, S., Ribeiro, A., Grannan, A., King, E. M. & Aurnou, J. M. 
2015 Laboratory-numerical models of rapidly rotating convection in planetary cores. Geo- 
phys. J. Int. 201 , 1-17. 

Davies, C.J., Gubbins, D. & Jimack, P.K. 2009 Convection in a rotating spherical fluid 
shell with an imposed laterally varying thermal boundary condition. J. Fluid Mech. 641 , 
335-358. 

Dowling, T. E. 1988 Rotating Rayleigh-Benard convection with fixed fluix boundaries. Sum¬ 
mer Study Program in Geophysical Fluid Dynamics, Woods Hole Oceanographic Institution 
Technical Report pp. 223-250. 

Greenspan, H. P. 1968 The Theory of Rotating Fluids. London: Cambridge University Press. 

Heard, W. B. & Veronis, G. 1971 Asymptotic treatment of the stability of a rotating layer 
of fluid with rigid boundaries. Geophys. Fluid Dyn. 2 ( 1 ), 299-316. 

Hide, R. 1964 The viscous boundary layer at the free surface of a rotating baroclinic fluid. 
Tellus 16 , 523-529. 

Hurle, D. T. J., Jakeman, E. & Pike, E. R. 1967 On the solution of the Benard problem 
with boundaries of finite conductivity. Proc. R. Soc. Lond. A 296 , 469-475. 

Johnston, H. & Doering, C. R. 2009 Comparison of turbulent thermal convection between 
conditions of constant temperature and constant flux. Phys. Rev. Lett. 102 (064501). 

Julien, K., Aurnou, J. M., Calkins, M. A., Knobloch, E., Marti, P., Stellmach, S. & 
Vasil, G. M. 2015 A nonlinear model for rotationally constrained convection with Ekman 
pumping. J. Fluid Mech., under review . 

Julien, K., Knobloch, E., Rubio, A. M. & Vasil, G. M. 2012a Heat transport in Low- 
Rossby-number Rayleigh-Benard Convection. Phys. Rev. Lett. 109 (254503). 

Julien, K., Rubio, A. M., Grooms, I. & Knobloch, E. 20126 Statistical and physical bal¬ 
ances in low Rossby number Rayleigh-Benard convection. Geophys. Astrophys. Fluid Dyn. 
106 (4-5), 392-428. 

Nayfeii, A. H. 2008 Perturbation methods. John Wiley & Sons. 

Nieves, D., Rubio, A. M. & Julien, K. 2014 Statistical classification of flow morphology in 
rapidly rotating Rayleigh-Benard convection. Phys. Fluids 26 (086602). 

Sprague, M., Julien, K., Knobloch, E. & Werne, J. 2006 Numerical simulation of an 
asymptotically reduced system for rotationally constrained convection. J. Fluid Mech. 551 , 
141-174. 

Stellmach, S., Lisciiper, M., Julien, K., Vasil, G., Cheng, J. S., Ribeiro, A., King, 
E. M. & Aurnou, J. M. 2014 Approaching the asymptotic regime of rapidly rotating 
convection: boundary layers versus interior dynamics. Phys. Rev. Lett. 113 (254501). 

Takeiiiro, S.-L, Isiiiwatari, M., Nakajima, K. & Hayashi, Y.-Y. 2002 Linear stability of 
thermal convection in rotating systems with fixed heat flux boundaries. Geophys. Astrophys. 
Fluid Dyn. 96 (6), 439-459. 


